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METHOD AND APPARATUS FOR COMPENSATION FOR 

NON-UNIFORM • ILLUMINATION 

packgrpunfl pf the lnv*n-tfen. 

This invention relates to compensation for non- 
5 uniformity in illumination of a surface in a system that 
takes information in optical form from the surface and 
converts such information to electronic form and in par- 
ticular to compensating for variations in illumination in 
systems for reading labels on articles, such as packages 
10 or letters , being moved at high speed along a conveyor 
belt past a video camera. 

In the shipping of parcels and packages, it is 
necessary, at various points in the process, to obtain 
various information that may be encoded on a label. Such 
15 information includes address information, shipper informa- 
tion, and the like. Encoded labels also have application 
in other fields, such as manufacturing. Uses for informa- 
tion encoded on a label may include, by way of example, 
sorting packages further along the conveyor belt, deter- 
20 mining the quantity of items going to a particular loca- 
tion, and recording the movement of packages so as to 
provide updates in response to customer inquiries. The 
rapid obtaining of information from a label may be facili- 
tated by an automated system for converting optically 
25 encoded information on the labels to electrically encoded 
information. 

In modern terminals for sorting of packages, 
packages are loaded on a conveyor belt moving at speeds of 
up to 100 inches per second and moved to various points 
30 for sorting. Such conveyor belts are ordinarily three to 
four feet vide. A common width is 42 inches. Packages 
may be loaded onto such conveyor belts anywhere across 
their widths, and in any orientation. Packages of differ- 
ing sizes may be shipped along a single conveyor belt. As 
35 a result when a package passes through a zone for reading 
of label information, the label may be at any point across 
the width of the conveyor belt and at any height within a 

range of heights above the surface of the conveyor belt. 

-1- 
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At. 

Objects of th e Invention. 

Accordingly, it is an object of this invention 
to provide a method and an apparatus for correcting sig- 
nals representative of information stored on a. surface in- 
5 optical form detected under conditions of nbri~unif ormity 
* of illumination of the surface.. r : . 

It is a further object of this invention to 
provide such a method .and apparatus for use in a system in 
which labels are placed on packages of varying heights. 
10 It is a further object of this invention to 

provide a method and apparatus for correcting signals 
representing label information in which labels are placed 
on packages that are moving on a conveyor belt. 

It is a particular object of this invention to 
15 provide such a method and apparatus for use where the 
information is stored on the surface in white areas and 
black areas. 

Additional objects and advantages of the inven- 
tion will become evident from the detailed description of 

20 a preferred embodiment below. 
Summary o f the Invention. 

A method for illumination compensation for use 
in a system for reading surfaces having at least black and 
white areas, includes the steps of storing reference 

25 reflected light intensity information corresponding to at 
least black and white areas in the form of the plurality 
of intensity values f each one of the intensity values 
having an address defined by a value of at least one 
variable, calculating from the reference reflected light 

3 0 intensity information more than one black threshold value 
and more than one white threshold value each of which 
threshold values is associated with a selected value of 
each of the variables. The method further includes the 
steps of receiving opera-tional reflected light intensity 

35 information in the form of more than one operational 

reflected light intensity value, each of the operational 
reflected light intensity values having an address defined 
by a value of each of the variables, remapping each of the 
operational reflected light intensity values to one of at 
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a preselected threshold minimum reflected light intensity 
value associated with the value of, each of the variables 
for the reflected light intensity value to a minimum level 
output and memory means for remapping any of thie reflected 
5 };icj:.t intensity values eouai to or greater than a prese- 
lected threshold maximum reflected light intensity value 
associated with the value of each of the variables for the 
reflected light intensity value to a maximum level output. 
Brief Description of the Drawings . 
10 Figure 1 is a simplified perspective view of a 

system for use with a method and apparatus according to 
the invention. 

Figure 2 is a block diagram of a circuit accord- 
ing to the invention. 
15 Figure 3 is a partial view of an calibration bar 

for use in connection with the invention. 

Figure 4 is a flow chart of the overall process 
of calibration in connection with a method and apparatus 
according to the invention. 

2 0 Figure 5 is a graph showing reference reflected 

light intensity data versus horizontal position before 
simplification. 

Figure 6 is a flow chart showing the steps of 
reflected light intensity data simplification. 
25 Figure 7 is a graph showing reference reflected 

light intensity data versus horizontal position after data 
simplification. 

Figure 8 is a graph showing output reflected 
light intensity values versus input reflected light inten- 

3 0 sity values as remapped according to the invention. 

Detailed Description of a Preferred Embodiment. 

With reference to Figure 1, there is shown an 
example of an apparatus for reading surfaces having infor- 
mation, such as reading label information on articles 
35 being moved on a conveyor belt. Such an apparatus may be 
used in connection with an illumination compensation 
apparatus of the invention. There is shown a section of a 
conveyor belt 60. Conveyor belt 60 is moving in the 
direction of arrow 61. Box 62 is on the conveyor belt 60. 

-5- 
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to discrete pixels. The reflected light intensity infor- 
mation is in the form of a series of discrete values or 
reflected light intensity values. Each discrete value 
corresponds to a particular position along a scan line. 
5 A second type, of information is information 

relating to the valvies.of the variables, or. variable value 
information, upon which the illumination depends. In the 
disclosed embodiment of the invention, there are two 
variables upon which the illumination depends. One such 

10 variable is height. Height information is transmitted 
along line 34 to a computer 36. 

The second such variable is horizontal position. 
The horizontal position of any item of reflected light 
intensity is preferably derived from a clock signal asso- 

15 ciated with the reflected light intensity information 

transmitted by the light-detecting device- A clock signal 
from the light-detecting device is transmitted along line 
3 8 to position logic unit 40. The light-detecting device 
also transmits a signal to the position logic unit 40 

20 representing the beginning of a new scan line. From the 
clock signal and the new scan line signal, position logic 
unit 4 0 provides a position signal. The unit may be an 
Altera 5032. The providing of such a signal is conven- 
tional in the art. The position signal is sent simulta- 

25 neously with each intensity value. The position signal 

counts from zero on each scan line across the illumination 
zone. The position signal value is changed a certain 
number of times in each scan line depending on detail of 
compensation in the apparatus. This is explained below in 

3 0 more detail where the look-up tables are discussed. 

computer 36 is shown in Figure 2. As noted 
above, an input of computer 36 receives height informa- 
tion. Computer 3 6 has the function of processing height 
information in the system of the invention. One output of 

35 computer 36 is coupled to an input of initial gain unit 

31. This output is a gain signal. In a preferred embodi- 
ment, this gain signal is an 8-bit digital signal. This 
gain signal is used as an adjustment for gain variation 
relating only to height. It has been found that reflected 

-7- 
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of two variables define the address of each reflected 
light intensity value. 

During calibration of an'-.apparatus according to 
the invention, the reflected light intensity value infor- 
5 mation represents light reflected from c surface having a 
known, or reference, optical property, ' Such reflected 
light intensity value information is referred to as refer- 
ence reflected light intensity value information. 

During operation of an apparatus according to 
10 the invention, the reflected light intensity value infor- 
mation represents light reflected t from a surface of un- 
known optical properties. Such reflected light intensity 
value information is referred to as operational reflected 
light intensity value information. 
15 The initial gain adjustment apparatus will now 

be explained in greater detail. An 8 -bit digital height 
signal is provided from the computer. The signal is input 
to a multiplier such as a CA 3338. The output of this 
logic unit is a gain-adjusted signal containing reflected 
20 light intensity value information. This is an analog 

signal. The amplitude of the output signal is determined 
by the height signal provided by the computer. In the 
preferred embodiment, a chip is provided that has 256 
possible gain adjustments. 
25 The height signal is generated by the computer 

depending on the height of the package. In one embodi- 
ment, a height sensor is employed, over a range of 28 
inches, with a sensitivity of one-half inch. Thus, in 
that embodiment 56 different height readings are possible. 
3 0 In another embodiment, the height sensor is used over a 
range of 3 6 inches, thus making 72 different height read- 
ings possible. The computer has a look-up table stored in 
memory of gain signal values corresponding to each possi- 
ble height value. The entries in the look-up table are 
3"5 obtained during calibration of the apparatus. The cali- 
bration of the apparatus is explained in detail below. 

The output of the gain adjustment unit is cou- 
pled to an analog-to-digital conversion unit. The unit 



-9- 



PCT/US92/07054 

changes 16 times in every sweep of the CCD across the 
width of the belt. The position signal is provided to at 
least two memory devices corresponding to left and right 
sides of the illumination zone. In a preferred embodi- 
5 ment, there are four channels provided. Two channels 

receive ref lected' licrht intensity value information corre- 
sponding to the left half of the illumination zone. Two 
other channels receive reflected light intensity value 
information corresponding to the right side of the illumi- 
10 nation zone. Thus the position value has 16 possible 
values , while 32 horizontal segments are provided. The 
same position signal corresponds to two position variable 
values, depending on which channel received the signal. 
The position value signal is thus a four-bit signal for 
15 each half of the illumination zone. Each time the posi- 
tion value changes, a new table is designated on each half 
of the illumination zone. 

In the preferred embodiment/ a CCD imager is em- 
ployed having 4,096 pixels in a scan line. Thus each 
20 table remaps 128 reflected light intensity values. 

The operation of the look-up tables in the 
memory devices is best understood with reference to Figure 
8. Figure 8 is a graph of output versus input for three 
horizontal segments at a single height range. Figure 8 
25 thus represents remapping of input to output for three of 
the 32 tables in a single height range. 

Line 802, which is a staircase curve, represents 
the remapping of the input values to the output values by 
one table. Line 802 corresponds to segment 2, which 
30 remaps values of pixels 128-255. All input digital inten- 
sity values less than a black level threshold value are 
remapped to a black output level. In the illustration of 
Figure 8, the black output level is zero. The black level 
threshold value for line 802 is approximately 102. All 
35 input values above a white threshold are remapped to a 
white output level. For the segments illustrated in 
Figure 8, the white output level is the maximum output 
level of the memory device. The white threshold for line 
802 is approximately 192. 

-11- 
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In such event, a reference reflected light intensity value 
would be obtained using a calibration bar having areas 
with whatever optical property represented the minimum re- 
flected light intensity value. This reference reflected 
5 light intiensity value would be designate a threshold 

minimum reflected light intensity value. The memory means 
would store a look-up table designed to remap any opera- 
tional reflected light intensity value equal to or less 
than the threshold minimum reflected light intensity value 
10 obtained from the reference to a minimum output level. 

Similarly, it is possible to use the invention 
in a context in which there is a maximum reflective light 
intensity value that corresponds to an optical property of 
a surface other than white. In such a system, there would 
15 be a surface used for calibration that had the optical 
property corresponding to the maximum reflected light 
intensity value used for determining a reference maximum 
reflective light intensity value. This reference reflect- 
ed light intensity value would be designated a threshold 
20 maximum reflected light intensity value. The memory 

device according to the invention would be programmed with 
a look-up table that remapped all operational reflected 
light intensity values equal to or greater than the thres- 
hold maximum reflected light intensity value to a maximum 
25 output level. 

It will be seen that in the general case, any 
operational reflected light intensity value falling be- 
tween the minimum reflected light intensity value thres- 
hold and the maximum light intensity value threshold would 
3 0 be remapped to an output value between a minimum output 

value and a maximum output value. For example, the values 
of operational reflected light intensity value falling 
between the minimum reflected light intensity value thres- 
hold and the maximum reflected light intensity value 
-35 threshold can be remapped, to an output value between a 

minimum output value and a maximum output value. However, 
the intermediate values could be remapped according to any 
desired function. 
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storage in a memory device. This is shown in box 516 in 
Figure 4. The last step is the step, shown in box 518 of 
Figure A, of storing the tables in the memory devices. 

In making these calibration adjustments , it is 
5 desirable to use a calibration bar 100 as shown in Figure 
3. Alternating white strips 102 aind black strips 104 are 
provided on the calibration bar 100 perpendicular to the 
length of the bar. For example, each strip may be 0.027 
inches in width. The calibration bar 100 is sufficiently 
10 long that it can rest on pegs 8 2 of Figure 1, and extends 
across the entire belt. It will be seen that for obtain- 
ing the reference reflected light intensity values, a 
surface or surfaces must be provided having both the 
optical property corresponding to the minimum reflected 
15 light intensity value and the optical property correspond- 
ing to the maximum reflected light intensity value. 

In the first step, shown as box 502, the attenu- 
ation gain is adjusted. The adjustments are made using an 
oscilloscope to monitor the levels. The goal in adjusting 
20 the attenuation gain is to assure that the black level is 
not below zero, and that the white levels are not saturat- 
ed. The calibration bar is placed at various selected 
heights. The step of adjusting the gain is repeated at 
each selected height. The gain value for each selected 
25 height is stored in the computer for later processing. In 
a preferred embodiment, the selected heights are 1 inch, 4 
inches, 8 inches, 12 inches, 16 inches, 20 inches, 24 
inches and 28 inches above the conveyor belt. In a system 
for reading surfaces having heights up to 36 inches, 
3 0 additional readings would be taken. 

The gain values set in the calibration process 
represent eight height points. In fact, in a preferred 
embodiment of the invention, there is provided height 
sensing means having the capability of detecting half-inch 
"35 height intervals. The points intermediate the eight 

selected heights are obtained by standard curve-fitting 
techniques. In a preferred embodiment, the identical 
curve-fitting technique and software that is used as de- 
scribed below in the step of simplification of the com- 

-15- 
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tially in the form of two lines. With reference to Figure 
5, the horizontal axis represents horizontal position 
across the illumination zone. The vertical axis represents 
reflected light intensity. In Figure 5, thirty-two data 
5 points , or pixels, have been chosen for display in th-.a 
figure. The intensity values shown have been measured by 
placing an input of an oscilloscope on an intensity output 
line from the light-detecting device. The intensity value 
numbers are arbitrary. 
10 The line designated by 602 in Figure 5 shows the 

white levels. As may be seen, there is substantial varia- 
tion in the white levels. In particular, the white levels 
drop on the edges of the illumination zone. The line 
designated by 604 represents the black levels. As may be 
15 seen, there is also considerable variation in the black 
levels. A profile similar to Figure 5 is obtained for 
each height interval at which the calibration bar is 
placed. The set of raw profiles are stored for further 
processing. Each raw profile is stored in a separate file 
20 in the computer. 

After collection and storing of the reference 
intensity values, there is an. optional step of checking 
the reference intensity value information data for evi- 
dence of hardware failure. A hardware failure in the 
25 system would make reading of labels with a reasonable 
degree of accuracy impossible. Such hardware failures 
include a failure of a portion of the illumination appara- 
tus , such as a bulb, or of a portion of the light sensing 
means, such as a failure of a CCD array. 
3 0 Checking the profile data for evidence of hard- 

ware failures is accomplished by computing a variance and 
standard deviation. The equation for the variance, var 
(X UM ), is the following: N n i 

v«r(x,...„) = E(*J-(*'Zv^ 

u * 

where X is the intensity value and N is the number of 
35 pixels in a profile. The standard deviation <r(X ( ..n) 
is defined by the equation of (Tfx,..^)- J v ^r(% t ..r\) 
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tion zone. The equation for finding the fitted curve from 



the raw profiles is the following: 



57* I** I* 3 * b, = Z*;y 




The source code for this is shown in CURVE-FIT 
C f which is included herewith in the appended software 
5 listings. 

Referring to Figure 6, there is shown a flow 
chart depicting the steps of the process of data simplifi- 
cation employing a preferred curve-fitting technique, in 
the step illustrated by the box labeled 702 on Figure 6, 
10 there is shown the step of allocating arrays. The alloca- 
tion of arrays is a conventional preliminary step to 
performing vector and matrix mathematics. The next step 
is the step of inserting the X and Y values in the matri- 
ces and arrays. The X values are the pixel numbers, which 
15 correspond to horizontal position in the illumination 

zone. The Y values are reflected light intensity values. 
This step is shown at box 704 in Figure 6. The next step, 
shown in box 706 in Figure 6, is the step of computing the 
X summations. The X summations are elements of the left- 
20 hand matrix in the above equation for the curve-fit. 

This step requires computation of summations of an array 
of X-values raised to various powers. A preferred soft- 
ware subroutine for this computation is shown under the 
file name "sum^a" in the appended source code listings. 
25 The next step in the process is the step, shown in box 708 
in Figure 6, of computing the XY summations. This step 
reguires taking two arrays of data, the array of X-values 
and the array of Y -values, and calculating a summation. 
The values in the array of Y-values are raised to increas- 
3 0 ing powers. A software subroutine for this step is shown 
under the file name n sum_ab" in the appended source code 
listings. The next step, shown at box 710, is transfer- 
ring the values obtained from the step of computing the X 
summations to a matrix. This is the matrix shown on the 
35 left hand side of the above equation for' the second-degree 
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After the step of processing the raw profile 
data into curve-fitted profile data, there are four sets 
of profile data. Each set of profile data represents a 
profile corresponding tD one of four ranges of height 
5 values in the illumination zoxixi. Each set of profile ha 
is next broken into 3 2 segments as shown by box 512 in 
Figure 4. The 32 segments represent the 32 horizontal 
segments of the illumination zone. Each of the 3 2 seg- 
ments corresponds to 128 pixels across the illumination 
10 zone. 

The simplified reference reflected light inten- 
sity values, or simplified profile data, as shown by box 
514 in . Figure 4, is then mapped to an inverse curve. In 
obtaining the inverse, each line is first rescaled by 
15 conventional techniques. 

The values contained in the inverse curve are 
used in programming the illumination compensation memory 
devices . The values must be placed in a format that can 
be accepted by the memory devices. In a preferred embodi- 

2 0 ment, an Intel MCS-B6 Hex file is created for each of the 

PR0M f s. The data stored is in the form of a table having 
256 elements. 

The data is then converted to a form which is 
compatible with the various channels into which the data 
25 is divided by the CCD units. For example r in a preferred 
embodiment, there are four channels. These four channels 
represent respectively odd and even pixels on the left and 
right sides of the belt. Accordingly, the data generated 
is divided into appropriate segments for each channel. 

3 0 Once the data is divided in this manner, the data is 

stored in computer memory in four parts. The data is 
transferred from the computer to the memory device. In a 
preferred embodiment, the memory device is an ultra-violet 
light erasable programmable read-only memory (PROM) . A 
35 Cypress 7C254 may be used, for example. In order to 

transfer data from the computer to such a device, a PROM 
burner is provided. The input of the PROM burner is 
coupled to a port of the computer so that data may be 
transferred from the computer to the PROM burner. In 

-21- 
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* File Name: moment. c . . 

* File Function; Routine to compute the average, standard deviation 

* and varaiance of the data passed. ; ' v \* ". 

* 

* Notes: Function was taken from Numerical Recipe in • C* see * 

chapter 13 pg 4 75 

* this function is a subset of the moment function. 

* History: 1988 Created 
/include <math.h> 

void moment ( data , m , n , index , ave , sdev , svar ) 
int m,n, index; 

f loar **da-ca, *ave , *sdev, *svar ; 

{ 

int j; 

float s; 

void nrerror () ; 

if (n <= 1) 
nrerror ( M n must be at least 2 in MOMENT") ; 

s=0 . 0 ; 

for (j = l; j <= n;j++) 
s data [index] [ j] ; 

*ave=s/n; 

*svar = 0.0; 

for(j = l; j <= n; 

*svar *=* fabs(s = data [ index] [j ] - (*ave)); 



*svar n; 
*sdev=sgrt (*svar) ; 
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curvefit (data, type, index, b) - .. - 

float * : *asta,*b; 

int index.; 

char tvpt:; 
{ ' - 

float *x, *y,**a,*xsum,*xysum,*col, **ai , d r power ; 
int i, j ,k,m,*indx; 
FILE *fp; 

x = vect or (J, VECMAX ) ; /* allocate array space */ 

y = vector (1 , VECMAX) ; ' 

xcuiD = vector (1, VECMAX ) ; 

xysum = vector ( 1, VECMAX) ; 

a - matrix(l,ARRAMAX, 1, ARRAMAX) ; 

ai = matrix (1, ARRAMAX, 1 , ARRAMAX) ; 

col = vector (1, VECMAX ) ; 

indx = i vector (1 , VECMAX) ; 

if (type — »p'j /* Dealing with profile gains */ 

a = VECMAX; ' 
else 

m * VECMAX/4; /* Dealing with SA gains */ 

if (type == «P») /* Profile Gains */ 

{ 

for(i = l; i <= n; i++) /* Get x values */ 
x[i] = i; 

■} 

else /* Sa gains */ 

{ 

forfi = l, j = l; i <= m; i++, j += 4) 
x[i] = j; /* Note Sa gains are increments 4 */ 

forfi = 1; i <= m; i++) /* Get y values */ 

y[i] = data [index] [i] ; 

power =1.0; /* Compute the x summations */ 

forfi = 1; i <= SUMAX; i++) 
xsum[i] - sum_a fx, m, power ++) ; 



power =1.0; /* compute the xy summations */ 

xysum[l] = sum^a (y,m, power ) ; 

forfi = 2; i < SUMAX; 
xysumfi] = sum_ab(x,y , m,power++) ; 
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* File Name: sum_a 
* 

* File. Function: _ To take, an; .array of data and find a ■j.uumation * 

raised to a variable power 

* Notes: 

* History: JSK 4-June-91 Created 
* 

***************************************^ 

float sum_a (a, limit, power) 
float *a , power; 
int limit; 

{ 

int i; 

float sum = 0.0; 

for(i = 1; i <= limit; i++j 
sum += pov(a[i] , power) ; 

return sum; 

} 

/***************** *********************************************** 

* File Name: sum ab 
* 

* File Function: To take two arrays of data and find a summation 

* raised to variable power 

* Notes: 
* 

* History: JSK 4-June-91 Created 
* 

***************************************************^ 

float sum_ab(a,'b, limit, power) 
float *a, *b f power; 
int limit; 

{ 

int i; 

float sum = 0,0; 

for(i - 1; i <= limit; i++) 
sum += (pow(a[i] , power) * b[i]); 

return sum; 
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/**************** 

* File lia^e: ludcmp.c 

* 

* File Function: Given an nxn matrix routine replaces it by the LU 
'* decomposition of a row vise per^a-cation of * 

itself. 

* Notes: Function was taken from Numerical Recipe in « C see * 

chapter 2 pg 4 3 

* History: 1988 Created 
* 

******************************* **** ************************* *** ** J 

/include <math.h> 

/define TINY 1.0e-20; 

void ludcmp(a, n, indx, d) 
int n,*indx; 
float **a,*d; 

{ 

int i , iaax , j , k ; 

float big , dum , sum , temp ; 

float * w , * vector ( ) ; 

void nrerror() , free_vector () ; 

w«vector(l,n) ; 
*d«1.0; 

for (i*=l;i<«n;i++) { 
big^O.o; 
for ( j<-n? 

if ( (temp=f abs (a[i] [ j] ) ) > big) big=temp; 
if (big — 0.0) nr error ("Singular matrix in routine 
LUDCMP"); w[i]=1.0/big; 
) 

for (j-l;j<-n;j++) { 

for (i«=l;i<j ;i++) { 
sum~a[i] [ j] ; 

for (k=l;k<i;k++) sum — a[i] [k]*a[k] [ j] ; 
a[i] [ j]=sum; 

} 

big=0.0; 

for (i=j ;i<=n;i++) { 
sum=a[i] [ j] ; 
for (k=l;k<j ;k++) 

sum — a[i][k]*a[k][j]; 
a[i] [ j]=sum; 

if ( (dum«w[i]*fabs(sum) ) >= big) { 
big=dua; 
imax=i ; 

} 

} 
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/* 

* File: genprom.c - - - 

* Auth: Robert K Pe)carek . 

* . Date: 03 Jan 1990 

* Copyright (c) 1990 United Parcel Service, 

* Purpose: Generate prom object files for SA correction tables. 

* Edit: Jim S Kunicki (JSK) added /define min marco to compile 

* under Unix 

* Changed progam to be a callable function 
*/ 

/define HALVES 2 
/define FIRST 0 
/define SECOND 1 
/define HZ 4 
/define PZ 16 
/define DATABITS 25 6 
/define SIZE 16 
/define VERSION 1 
/define REVISION 0 

/define min (a, b) (((a) < (b) ) ? (a) : (b) ) 



/include <stdlib.h> 
/include <stdio.h> 
/include <string.h> 
/include <ctype.h> 
/include <math.h> 
/ifdef MSDOS 
/include <process.h> 
/endif 

extern void main (int argc, char **argv) ; 

static FILE *GetInputFile(int argCount, char **stringVector) ; 
static void GetFSOpsRND (FILE *fp, int *full, float +preScale, int 
♦quantize, 

int *inputScale) ; 
static void BuildRawoutOff set (FILE *fp, float *offl, float *ravi, 

float *off2, 

float *rav2) ; 

static void CalcProm (float *offl, float *ravl , 

int full, float preScale, int quantize, 

int inputScale, 

unsigned char ^result, float *scale) ; 
static void WritelntelHex (unsigned char *result, int aLen, char 
*outFileName) ; 

static void WriteText (unsigned char *result, char *outFileName) ; 
static void WriteScaleFactors (float *scale, int aLen, char 
♦outFileName) ; 

******************** A*******************************************, 
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/* Write prom data to files */ 

WritelntelHex (Sresults [FIRST] [0] [0] [0] , £Z*HZ*DATABITs 
"AB.FRM") ; . ' 

WritelntelHex (&re suits [SECOND] [0] [0] [0] , P2*HZ*DATABITS 
"CD.PRM") ; ' 

/* Write prom data to files */ 

WriteText ( ^results [ FIRST ] [ 0 ] [ 0 ] [ 0 ] , " AB . TXT" ) ; 
WriteText (firesults [ SECOND] [0] [0] [0] , "CD. TXT" J ; 

/* Write out scale factors */ 

WriteScaleFactors (& scale [FIRST] [0] [0] , PZ*HZ*HALVES 
"SF.TXT") ; . ' 

/* Acknowledge completion */ 

fprintf (stderr, "Done generating prom files and scale 
factors \n") ; 

} 

/* 

* GetlnputFile - verify argument count and retrieve file name 
V 

/* */ 

static FILE *GetInputFile() 

{ 

FILE *tf ; 

/* Return invalid argument count */ 
/********************** 

if (argCount i= 2) 

{ 

return NULL; 

} 

********** ****************** ^^^^^t^^^^^^^^^^^A^^^^^^^^^^^^^^^^^^^ 

/* open file note hard coded for now */ 
if ((tf = f open ("DATA. OUT" , "rt") ) = NULL) 
{ 

perror ("") ; 
return tf; 

} 

/* 

* GetFSOpsRND - read scalers from input file in designed format 

/* */ 

/ifdef SPARC 

/define _MAX_PATH 144 /* Note the was a MICROSOFT "C" define 

within */ 

#endif /* stdlib.h, SUN "C" is not supported */ 
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perror ) ; 

exit(l) ; - - ^ 

• } 

} 

/* 

* BuildRavOutOf f set - read raw outputs and offsets from input file 
*/ 

/* */ 

static void BuildRavOutOf f set (f p, of f 1 , ravi , of f 2 , rav2 ) 
FILE *fp; 

float *offl,*ravl, *of f2 , *rav2; 

< 

/define NUMPARMS 2 

/define LONGBUFFER 300 

/define STARTPOS 23 

static char MatchString[ ] = "Diag" ; 

char buff er [LONGBUFFER] , *rc, *bufferPtr; 

float *ptrl, *ptr2; 

int hz f pz, count, buf f erof f set; 

/* Find first line containing the identifying string */ 
do 

{ 

rc = fgets (buffer, LONGBUFFER , fp) ; 

} while ((rc i= NULL) strncmp (buf f er , MatchString, 

stxlen (MatchString) ) ) ; 

if ( r c = NULL) 

{ 

perror ("No data found"); 
exit(l) ; 

} 

/* For all height zones in the file, and not at end of file 

*/ 

for (hz = 0; (hz < HZ) £& (rc != NULL); ++hz) 

{ 

/* For both sets of input data */ 

for (ptrl = offl, ptr2 = off 2 + PZ - 1, count - 0; 
count < NUMPARMS; 

ptrl = ravi, ptr2 = rav2 + pz - l, ++count) 

{ 

/* Locate beginning of data on line */ 
bufferPtr - ^buffer [STARTPOS] ; 

/* Get the data for the left side */ 
for (pz = 0; pz < PZ; ++pz, ++ptrl) 
{ 

sscanf (buff erPtr , " %f,%n", ptrl, 
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if 



(tempStore 



> 



(float) full) 



tempStore 



(float) full; 



tempStore 



(float) floor (tempStore / 

(float) quantize + (float) 0.5) \ * 

(unsigned char) tempStore ; 



(float) quantize ; 



*result++ 



} 



/* Increment pointers */ 

ravl++; 

offl++; 



} 



} 



/* 



* WritelntelHex - put character data out to given file .in the 
intel hex 



static void WritelntelHex (result , aLen, outFileName) 
unsigned char ^result; 
int aLen; 

char *outFileName ; 

{ 

FILE *fp; 

int count, tempLen, chkSum, dataCount; 
/* open file */ 

fp «= f open (outFileName, "vt") ; 

/* For each data group of SIZE */ 
count - 0; 
while (aLen > 0) 

{ 



/* Put out intel Hex line */ 
tempLen = min(SIZE, aLen); 

chkSum = 0x7 000 - tempLen - (count / DATABITS) - (count 



k DATABITS) ; 

fprintf(fp, ":%02X%04X00", tempLen, count); 



/* Put out data */ 

for (dataCount = 0; dataCount < tempLen; ++dataCount) 
i 

fprintf(fp, "%02X", *result) ; 
chkSum — *result++; 

} 

/* Put out checksum */ 

fprintf(fp, "%02X\n", chkSum % DATABITS); 

aLen -= tempLen; 
count tempLen; 



format 



*/ 



/* 
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fprintf(fp, "Channels A/B\n") ; '" 

/* Write all ccale factors */ 

for (count ~- 0: count < aLen/2; ++count) 

{ 

fprintffxp, "%2i *2i *4.1f\n" , count / PZ, count % pz 

*scale++) ; 

} 

fprintf(fp, "Channels C/D\n") ; 

for (count = 0; count < aLen/2; ++count) 

{ 

fprintf(fp, "*2i *2i *4.if\n", count / pz, count % pz 
*scale++) ; ' 

} 

f close (fp) ; 

} 
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tween the white output level and the black output level, 
multiplied by the quotient of said operational light 
intensity value less the black level threshold divided by 
the absolute value of the difference .between the black 
0 level threshold and the white level threshold. 

3. A method as recited in claim 1, wherein - ; 
said black output level is zero, and said white output 
level is a maximum output level. 

4. A method as recited in claim 1, wherein 
said step of storing comprises storing reference reflected 
light intensity information corresponding to at least 
black and white areas in a plurality of profiles, each one 

5 of said profiles displaying said reference reflected light 
intensity information plotted against one of said vari- 
ables, each other of said variables being a selected 
constant for each profile. 

5. A method for illumination compensation in a 
system for processing of signals representing reflected 
light intensities from a surface, comprising the steps of: 

(a) receiving reflected light intensity 

5 information in the form of a plurality of intensity val- 
ues ; 

(b) receiving, associated with each one of 
said reflected light intensity values, a value of at least 
•one variable; 

(c) remapping each of said reflected light 
intensity values to one of a plurality of compensated 
reflected light intensity values, including the steps of 
remapping any of said reflected light intensity values 
equal to or less than a preselected threshold minimum 
reflected light intensity value associated with the value 
of each of said variables for said reflected light inten- 
sity value to a minimum level output and remapping any of 
said reflected light intensity values equal to or greater 
than a preselected threshold maximum reflected light 
intensity value, associated with the value of each of said 
variables for each reflected light intensity value, to a 
maximum level output. 
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20 threshold maximum reflected light intensity value associ- 
ated with the value of each of said variables for each 
reflected light intensity value,. -.to a maximum level out- 
put. 

9 . An apparatus as recited in claim 8 , further 
comprising mer -ory means for remapping any of said ref'ecv.- 
ed light intensity values less than said preselected 
threshold maximum light intensity value associated with 
5 the value of said variables for said reflected light 

intensity value and greater than said preselected thres- 
hold minimum reflected light intensity value corresponding 
to a value of each of said variables for said reflected 
light intensity value, linearly to a value greater than 
10 said minimum level output and less than said maximum level 
output . 
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